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ABSTRACT 



Aims. A number of recent works have suggested that the period-luminosity (PL) relation for the Large Magellanic Cloud (LMC) 
Cepheids exhibits a controversial nonlinear feature with a break period at 10 days. Therefore, the aim of this Research Note is to test 
the linearity/nonlinearity of the PL relations for the LMC Cepheids in BVI^JHR, band, as well as in the Wesenheit functions. 
Methods. We show that simply comparing the long and short period slopes, together with their associated standard deviations, leads 
to a strictly larger error rate than applying rigorous statistical tests such as the F-test. We applied various statistical tests to the current 
published LMC Cepheid data. These statistical tests include the F-test, the testimator test, and the Schwarz information criterion 
(SIC) method. 

Results. The results from these statistical tests strongly suggest that the LMC PL relation is nonlinear in BVIcJH band but linear in 
the Ks band and in the Wesenheit functions. Using the properties of period-color relations at maximum light and multi-phase relations, 
we believe that the nonlinear PL relation is not caused by extinction errors. 
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' Recently, Fouque et al. (I2007I I have derived the Galactic Cepheid 
period-luminosity (PL) relation with several different tech- 
niques, including parallax measurements (from Hipparcos and 
HST), variants of the Baade-Wesselink method, and distances 
inferred from open clusters. We point out that such an ap- 
proach has been applied before in Ngeow & Kanbur (I2004I I and 
Groenewegen et al. ( |2004| l. In addition, Fouque et al. ( I2007I I 
also derive Large Magellanic Cloud (LMC) PL relations in the 
BVRJcJHKs band, and refer to the work of Sandage et al. 
( 1200 4). which suggests a possible change of slope for the LMC 
PL relation at 10 days. In fact, there are several other papers on 
the topic of nonlineaiQLMC PL relations (see Kanbur & Ngeow 
120041 [loos', Kanbur et al. 2007j Ngeow et al. 12005 ; Ngeow & 
Kanburg006a 2006b, Koen et al. l^OTl l. 

These previous works concentrate on the Vic band (Kanbur 
& Ngeow l2004l 120061 Ng eow & Kanbur'2006b) or V band only 
(Ngeow & Kanbur l2006al Kanbur et al. 2007), with data mosdy 
from the OGLE (Optical Gravitational Lensing Survey, Udalski 
et al. |I9991 l database. For the JHKs band PL relations, Ngeow et 
al. (I2005I I investigated possible nonlinearities using the 2MASS 
data from Nikolaev et al. ( 120041 ) that cross-correlated with the 



MACHO LMC Cepheids, and a random-phase coiTection to de- 
rive the mean magnitudes of these 2MASS data. 

Our motivation for this Research Note is to extend the previ- 
ous work in BVIcJHK, band, using the LMC Cepheid data from 
Fouque et al. (2007), with various rigorous statistical tests. The 
JHK, band data used in Fouque et al. (2007 ) are the 2MASS data 
matched to the OGLE Cepheids, and the mean magnitudes are 
derived using the method presented in Soszyiiski et al. (2005), 
which is dilferent from the data used in Ngeow et al. (2005). It 
is important to test the nonlinearity results in JHK, band results 
with different Cepheid samples and different methods deriving 
the JHKs mean magnitudes. 



Send offprint requests to: C. Ngeow 

' By nonlinearity we mean that the PL relation can be broken into 
two relations, with a break period adopted at 10 days. 



As emphasized in Ngeow & Kanbur ( |2006a! ). statistical tests 
are needed to test and detect the existence of the nonlinear PL 
relation. We also point out that in searching for nonlinearity or 
a change of slope at 10 days, the method of comparing the short 
and long period slope with their associated standard deviations 
is more prone to error than applying a statistical test, such as the 
F-test as indicated by the following, purely analytical example. 
A statement such as the "the slope is x + dx" means that the 
probability that the slope is in the interval (x - 6x, x + 6x) is 
1 - a, where a is the desired significance level. Then if A is the 
event that the calculated short period slope is wrong and B is 
the event that the calculated long period slope is wrong, we have 
F(A) = a and P{B) = a. Then in comparing the short and long 
period slopes using just their calculated standard deviations, the 
probability of at least one mistake is P(A U B) - P(A) + P(B) - 
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P(A n B) ^ 2a - a^. If I > a > Q, then 2a - a'^ > a. If the 
F-test or any other statistical test is carried out to the level of 
significance a, then this states that the probability of making an 
error in just comparing long and short period slopes through their 
standard deviations is greater than the probability that the F-test 
makes a mistake. In essence the F-test compares both short and 
long period slopes (as well as the zero-points) of the nonlinear 
PL relations simultaneously. 

2. Data & Results from Statistical Tests 

The BVIcJHKs band LMC Cepheid data were kindly provided 
by P. Fouque. It is exactly the same dataset used in Fouque et al. 
( 12007) . We do not include the Rc band data because the num- 
ber of Cepheids in Rc band data is much smaller than in other 
band and some of them have a small number of data points per 
light curve. In addition to the BVIcJHKs band, we also include 
the two Wesenheit functions considered in Fouque et al. (200^, 
namely Wu and Wv,. We fit a linear PL relation to these data and 
obtain identical PL relations as presented in table 8 of Fouque 
et al. (120071) . The statistical tests we employed in this study in- 
clude the F-test (Kanbur & Ngeow 2004; Ngeow et al. 120051 ). 
the testimator test, and the Schwarz information criterion (SIC) 
method (Kanbur et al. 2007). Details regarding the formalism 
and description for these statistical tests are given in the above 
references and will not be repeated here. 

To illustrate the difficulty of visualizing the nonlinear PL re- 
lation, if it truly exists, and the need for rigorous statistical tests 
to detect such nonlinearity, we simulate the J band PL relation 
using the method detailed in Ngeow & Kanbur ( |2006al l. The in- 
put linear J band PL relation is taken from Fouque et al. (2007), 
and the input nonlinear PL relation in the simulation is adopted 
from Table [T] Figure [T] compares the J band PL relation from 
the real data and the two simulations. The three PL relations in 
Figure[T]look similar and linear by eye. However the PL relation 
in the bottom panel is constructed from a nonlinear PL relation. 
As pointed out in Ngeow & Kanbur (2006a), the difficulty of vi- 
sualizing such a nonlinear PL relation is due to the existence of 
intrinsic dispersion of the PL relation caused by the finite width 
of the instability strip. 

2.1. The F -Test Results 

For the F-test, the null hypothesis (Ho) is a single regression 
line is sufficient, while the alternate hypothesis (Ha) is that two 
regression lines separated at 10 days are needed to fit the data. 
In Table[T] we present the results from the F-test, which include 
the fitted PL relations for the short (log P < l.O) and long period 
Cepheids, the F values and the probability, piF), under the null 
hypothesis. As in our previous work, the threshold for p{F) was 
set to be 0.05 (corresponds to 95% confident level). For a large 
sample (A^ > 100), F ~ 3 at p(F) = 0.05. Hence our F-test 
results indicate that the LMC PL relation is not linear in BVIcJH 
band but linear in KsWbiWyi band. Note that some of the slopes 
and/or zero-points in Table [T] appear to be consistent between 
the short and long period Cepheids, but this does not negate the 
nonlinearity of the PL relation as pointed out in the Introduction. 

2.2. The Testimator-Test Results 

Briefly, the testimator test requires the Cepheid sample to be di- 
vided into a number of sub-samples, after the sample has been 
sorted according to the periods. The slope of each sub-sample 
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Fig. 1. The top panel shows the J band data from Fouque et al. 
(2007). Both of the middle and lower panels are from simulation. 
The simulation in the middle panel uses an intrinsic linear PL 
relation as input, while the PL relation in the bottom panel is 
simulated from an intrinsic nonlinear PL relation. The F-test as 
described in Section 2.1 returns F = 0.68 and F = 7.62 for the 
PL relations in the middle and bottom panels, respectively. The 
F-test coiTectly identified the PL relations that are intrinsically 
linear and nonlinear, respectively. 



is compared to the slope from the previous sub-sample (except 
the first sub-sample) using the t statistical test. The null hypoth- 
esis for the f-test is that the two slopes are statistically equal, and 
the alternate hypothesis is that they are not. If the null hypoth- 
esis is rejected, a new slope, called the testimator, is calculated. 
The ratio, k, of the observed f values and the critical t value is 
calculated for each sub-sample (see Kanbur et al. l2007l for more 
details). The case of k > I indicates that the null hypothesis can 
be rejected for a given sub-sample, and the slope for this sub- 
sample is statistically different to the slope in the previous sub- 
sample. This ensures that the testimator tests only the slope and 
not differences in zero point. In Figure |2] we present the k val- 
ues of each of the sub-samples under the testimator test. From 
the figure, it can be seem that in the BVIcJH band, the slopes 
change (with A: > 1) for sub-samples that bracket the assumed 
break period at 10 days and/or the sub-samples with longer pe- 
riod Cepheids. In contrast, the slopes do not change statistically 
in the sub-samples for band. 

2.3. The SIC-Test Results 

For the SIC test, the null hypothesis is taken to be a linear re- 
gression model, while the alternate hypothesis is a nonlinear re- 
gression model with a break period at Pq. This break period is 
varied over the entire period range and likelihoods under the null 
and alternate hypotheses are calculated. We look for values of 
P() for which the likelihood under the alternative hypothesis is 
greater than that under the null hypothesis. The model with the 
lowest S IC value is the preferred model. Figure [3] summarizes 
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Table 1. F-test Results of the LMC PL Relations. 
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The subscripts s and l are for the short (log P 
dispersion of the fitted PL relations. 



1.0) and long period Cepheids, respectively, while a, b and cr are the slope, zero-point and 



1 1.5 

log (P) 



1 1.5 

log (P) 



Fig. 2. The results from the testimator test. The horizontal bars are the k values for each sub-samples. The size of the bars indicates 
the period range covered in each sub-samples. The number of data points in each sub-samples ranges from 45 to 88. The dashed 
lines are for the case k= I. 
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Fig. 3. The results from the SIC test. The thick horizontal lines (since it is independent of log Pq) and the "curves" are for the null 
and altemate hypothesis, respectively. The vertical dashed lines indicate the adopted break period at 10 days. 
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the results from this SIC test. In BVIJH band, the SIC test finds 
evidence that there is a range of Pq where the ahernate hypoth- 
esis is a preferred model. This range of Pq includes the adopted 
break period at 10 days (except in the / band, however the range 
of the break period is still close to 10 days). The existence of a 
range of Pq and the difficulty of pin-pointing the break period 
in the SIC test is mainly due to the finite width of the instability 
strip (see more detailed discussion in Kanbur et al. |2007| l. The Ks 
band results do not show any preferred alternate models around 
10 days. The SIC results for the two Wesenheit functions also do 
not prefer the alternate hypotheses across all Pq. Interestingly, 
the HKs band imply that at log Pq ~ 0.5 the alternate hypothesis 
show a smaller value of SIC than the null hypothesis. Currently 
there is no explanation to account for this. 

2.4. Tests for the JHK, Band PL Relations with Additional 
Data 

In contrast to other band, the JHKs data from Fouque et al. 
( 1200 7) only consist of the OGLE Cepheids. Therefore, we in- 
clude additional JHKs band Cepheid data from Persson et al. 
( I2004I I to the sample, as suggested by the referee. As in Kanbur 
& Ngeow (I2OO6I 1 and Fouque et al. (I2007I I. we only include 
Cepheids with log(f ) < 1.8 in the sample. Further we remove 
Cepheid HV 12765 from the sample as suggested in Persson et 
al. C2004 ). Since the JHKs data from F ouque et al. (200?) is in 
2MASS system and the Persson et al. (I2004II JHK ^ data is in 
LCO system, we convert the Persson et al. (120041 ) JHKs data 
to the 2MASS system using a shift of -0.02mag. as stated in 
Fouque et al. (120071 ). After applying the extinction correction, 
the linear JHKs PL relations with the combined 615 Cepheids 
are: 

m{J) = -3.121(+0.018)log(f)-H 16.263(+0.014), 
m{H) = -3.228(+0.016)log(P)+ 16.047(+0.013), 
m(Ks) = -3.249(+0.018)log(P)-t- 16.001(±0.015), 

with dispersion of 0.130, 0.1 17 and 0.133 respectively. 

The F-test results for the combined JHKs data is presented 
in Table |2] while the results from the testimator test and the SIC 
test are collectively summarized in Figure |4] All three statistical 
tests again found strong evidence of nonlinear JH band PL rela- 
tions, at the assumed break period of 10 days, for the combined 
Cepheid data. For the Ks band, the F-test shows that the PL re- 
lation is marginally linear, with the support from the testimator 
and SIC tests that the nonlinearity is not detected at the break 
period of 10 days. 

3. Conclusion & Discussion 

Combining the results from the three statistical tests presented 
in the previous sections, we find that there is strong statistical 
evidence to suggest the LMC PL relation is nonlinear in the 
BVIcJH band but linear in the KsWh/Wri band. Including ad- 
ditional data from Persson et al. (2004 ) for the JHKs band does 
not alter the results as well. We have to emphasize that both of 
the testimator and SIC methods are applied to the BIcJHKs band 
and the Wesenheit functions for the first time, in contrast to the 
V band data that has been studied in Kanbur et al. (120071 ). 

The nonlinear LMC PL relation has been found from a 
Cepheid sample that consists of OGLE Cepheids only (Kanbur 
& Ngeow 120041 ). To extend the OGLE sample, mostly at the 



long period end, Ngeow & Kanbur ( |2006al ) included various ad- 
ditional data from literature (see table 1 of Ngeow & Kanbur 
2006a), and again found strong evidence of nonlinearity of 
the LMC PL relation. In this Research Note, results using the 
Fouque et al. (120071 ) data alone, and with additional data from 
Persson et al. (120041 ). further supports the conclusion given in 
Ngeow & Kanbur (I2006al ): that the sample selection does not 
play an important role in detecting the nonlinear LMC PL rela- 
tion. However, Fouque et al. (120071 ) have suggested that the mix- 
ture of data used in previous work may lead to the nonlinearity 
seen in the statistical tests. This may certainly be the case and 
the analysis of a homogeneous sample, such as that provided by 
the "LMC shallow survey" (Fouque 2007; Gieren 2007 - private 
communication) is desirable. 

The nonlinearity of the PL relation that is seen in the opti- 
cal and JH band but not in the reddening insensitive Ks band 
and the Wesenheit function may suggest that extinction is the 
cause of the nonlinearity. However, extinction is not the only 
explanation and there is some evidence against the hypothesis 
of extinction errors as a cause for the apparent nonlinearity. The 
linearity of the Ks band PL relation, as compared to other shorter 
wavelength PL relations, is expected from black-body arguments 
(Ngeow & Kanbur I2006al ). Simply speaking, the temperature 
variation dominates the luminosity variation in the optical, and 
extends to JH band for Cepheid-like temperatures. But in the Ks 
band the luminosity variation is dominated by the radius vari- 
ation of Cepheid variables. The proposed mechanism that may 
cause the nonlinear PL relation, the interaction between the hy- 
drogen ionization front and the stellar photosphere (Kanbur & 
Ngeow 2006 ), will only affect the temperature variation and not 
the radius variation. The linearity of the Wesenheit functions is 
also not a surprise, and has been studied and discussed in Ngeow 
& Kanbur (i2005i ) and in Koen et al. (120071 ). and will not be re- 
peated here. 

Since the additional data used is mainly at the long period 
end, the possibility remains of systematic errors in reddening as 
a function of period. However we note that a reddening error as 
a function of long period LMC Cepheids would also change the 
observed properties of LMC Cepheids at other phases. A redden- 
ing error as a function of period such that LMC Cepheids obey a 
linear PL relation at mean light would force the LMC Cepheids 
to have a period-color relation such that they get bluer or hot- 
ter at maximum light as the period increases (see Figure |5] for 
a schematic illustration). This is in stark contrast to the behav- 
ior of Galactic Cepheids and long period LMC Cepheids, which 
are known to have a flat period-color relation at maximum light 
(Code [19471 Simo n et al [19931 Kanbur & Ngeow [20041 [2056l 
Kanbur et al'2004'). Moreover, it is difficult to explain, theoret- 
ically, how a Cepheid could get hotter at maximum light as the 
period increases. 

The PL relation at phase ~ 0.8 described in Ngeow & Kanbur 
(|2006bl ) presents clearly the dramatic nature of the nonlinearity 
at 10 days and the dynamic nature of the PL relation as a func- 
tion of phase. It is difficult to reconcile this behavior as being 
due to sampling errors and/or reddening errors. It is worth to 
point out that the mean light PL relation used in the literature is 
an average of the PL relations in all phases (Kanbur & Ngeow 
2004; Kanbur et al. 120041 Ngeow & Kanbur 2006b ). nonlineai-- 
ity of the PL relations at certain phases will certainly affect the 
linearity/nonlinearity of the mean light PL relation. 

We have to remind that the data used in this study (and 
in most of our previous work) were published data that have 
been corrected for extinction using the "state-of-the-art" and 
well-developed methodology. If extinction error is believed to 
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Table 2. F-test Results of the JHK, PL Relations with additional data from Persson et al. (! 2004| l. 
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Fig. 4. Test results for the testimator test (left panel) and the SIC test (right panel) for the JHKs band PL relations with additional 
data from Persson et al. (i2004) . See Figure|2]and[3]for the meaning of the lines and curves. 
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Fig. 5. Schematic illustration for the argument with PC(max) re- 
lation. Top panels show the observed PC relations (after cor- 
rected for extinction) at mean (top-left panel) and maximum 
(top-right panel) light. Bottom panels show that if additional ex- 
tinction as function of period to make the mean light PC relation 
linear, then the same extinction will cause the colors at maxi- 
mum light get bluer as period increases, which are against ob- 
servation and theoretical expectation. 



be the cause of nonlinear PL relations, it would imply that the 



extinction correction done previously in the literature is incor- 
rect and/or incomplete. This would affect the previous work that 
using these extinction corrections, and those results need to be 
revised in future work. 
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